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ABSTRACT 


Columns  are  the  key  load-bearing  elements  in  frame  structures  and  exterior  columns  are  probably  the 
most  vulnerable  structural  components  to  terrorist  attacks.  Column  failure  is  normally  the  primary  cause 
of  progressive  failure  in  frame  structures.  A  high-fidelity  physics-based  computer  program,  LS-DYNA  was 
utilized  in  this  study  to  provide  numerical  simulations  of  the  dynamic  responses  and  residual  axial 
strength  of  reinforced  concrete  columns  subjected  to  short  standoff  blast  conditions.  The  finite  element 
(FE)  model  is  discussed  in  detail  and  verified  through  correlated  experimental  studies.  An  extensive 
parametric  study  was  carried  out  on  a  series  of  12  columns  to  investigate  the  effects  of  transverse 
reinforcement  ratio,  axial  load  ratio,  longitudinal  reinforcement  ratio,  and  column  aspect  ratio.  These 
various  parameters  were  incorporated  into  a  proposed  formula,  capable  of  estimating  the  residual  axial 
capacity  ratio  based  on  the  mid-height  displacement  to  height  ratios. 


1.  Introduction 

Columns  are  the  key  load-bearing  elements  in  frame  structures. 
Exterior  columns  are  probably  the  most  vulnerable  structural 
components  to  terrorist  attacks.  Column  failure  is  normally  the 
primary  cause  of  progressive  failure  in  frame  structures.  However, 
current  knowledge  in  the  evaluation  of  residual  capacity  of  a  blast 
damaged  reinforced  concrete  column  remains  limited.  A  better 
understanding  of  residual  capacity  in  columns  would  aid  in  the 
prediction  of  the  overall  performance  of  buildings,  its  resistance  to 
progressive  collapse  and  determining  the  stability  of  damaged 
buildings  especially  during  search  and  rescue  operations. 

Single-degree-of-freedom  (SDOF)  analysis  from  blast-resistant 
design  guidelines  [1-3]  provides  engineers  with  simplified 
analytical  methods  to  assess  blast  damage  of  RC  columns.  Although 
these  simplified  methods  are  quite  useful,  three-dimensional 
analysis,  in  contrast,  provides  a  more  in-depth  understanding  by 
incorporating  all  aspects  of  the  response  of  concrete  structures 
subjected  to  blast  effects. 

A  study  by  Hayes  et  al  [4]  suggests  that  the  proper  application  of 
current-practice  seismic  detailing  for  high-seismicity  regions  can 
reduce  vulnerability  to  blasts  and  progressive  collapse.  One  of  the 
aims  of  this  study  is  to  quantify  this  improvement. 

A  three-dimensional  nonlinear  FE  analysis  utilizing  the  LS-DYNA 
software  [5]  is  performed  for  the  numerical  simulations  of 
this  research.  The  FE  model  is  validated  through  correlated 
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experimental  studies.  The  validated  FE  model  was  then  analyzed 
under  simulated  blast  loads  and  investigations  were  carried  out  on 
the  dynamic  responses  and  residual  axial  capacities  of  the  columns. 
An  extensive  parametric  study  was  carried  out  on  a  series  of  12 
columns  to  investigate  the  effect  of  the  transverse  reinforcement 
ratio,  long-term  axial  load  ratio,  longitudinal  reinforcement  ratio, 
and  column  aspect  ratio  on  the  column  responses. 

2.  Finite  element  model 

The  explicit  nonlinear  FEM  program  LS-DYNA  [5]  was  utilized  in 
this  study  because  of  its  proven  effectiveness  in  geometric 
modelling  and  impact  analysis.  The  description  of  modelling 
includes  blast  loadings,  the  structural  geometry,  relevant  material 
models,  application  of  loads  and  analysis  procedures. 

2.1.  Blast  loadings 

An  exterior  explosion  to  the  building  generates  four  types  of 
loads  as  shown  in  Fig.  1 :  impact  of  primary  fragments,  impact  of 
secondary  fragments,  overpressure,  and  reflected  pressure.  The 
study  reported  within  this  paper  is  restricted  to  the  effects  of 
overpressure  and  reflected  pressure  on  the  target  from  an  explo¬ 
sion.  As  the  overpressure  wave  strikes  on  the  front  face  of  a  closed 
target,  reflected  pressure  is  instantly  developed,  and  this  is  the 
most  destructive  aspect  of  blast  loading  on  a  structure.  In  this  study, 
the  explosion  centre  is  assumed  at  the  mid-height  of  a  column, 
while  the  surface  is  assumed  to  be  the  reflected  surface.  The  loading 
at  different  points  on  the  front  surface  of  the  column  for  a  given 
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Fig.  1.  Blast  loadings  on  the  first  fioor  column  during  a  close-in  explosion  [6]. 


charge  and  standoff  distance  is  computed  by  LS-DYNA  [5]  with  the 
built-in  CONWEP  blast  model,  which  relates  the  reflected  over¬ 
pressure  to  the  scaled  distance  and  also  accounts  for  the  angle  of 
incidence  of  the  blast  wave  [7].  Blast  incidents  in  recent  years 
showed  that  most  of  the  terrorist  attacks  on  public  structures  were 
explosions  within  short  standoff  distance  (<10m).  Thus  in  this 
study,  the  standoff  distance  is  assumed  to  be  5  m.  Considering  the 
limitation  of  the  weight  of  the  explosive  which  can  be  obtained  in 
a  particular  region,  a  maximum  weight  of  1  ton  equivalent  TNT  is 
presumed.  Various  charges  weighing  between  0  and  1  ton  are  used 
to  obtain  different  levels  of  damage  to  the  columns. 

2.2.  Structural  geometry  modelling 

Fig.  2  shows  the  three-dimensional  model  of  the  column.  Eight- 
node  solid  hexahedron  elements  are  used  to  represent  concrete. 
The  reinforcing  bars  are  modelled  explicitly  using  two-node 
Hughes-Liu  beam  elements  connected  to  the  concrete  mesh  nodes. 
The  nodes  that  link  the  concrete  and  reinforcement  mesh  are 
shared  and  these  shared  nodes  are  not  able  to  slip.  As  a  result,  they 
form  a  perfect  bond  due  to  this  assumption  of  complete  compati¬ 
bility  of  strains  between  concrete  and  steel.  The  restraint  at  the 
upper  end  of  the  column  provided  by  secondary  floor  beams  and 
slabs  is  modelled  as  a  stiff  block,  while  the  bottom  end  restraint  is 
modelled  as  a  flxed  support.  A  rigid  plate  which  is  only  allowed  to 
move  in  the  vertical  direction  is  attached  to  the  top  end. 


Concrete  mesh  Rebars  inside  concrete 


2.3.  Material  models 
23 A.  Concrete 

The  flnite  element  code  LS-DYNA,  which  is  used  in  this  research, 
contains  several  material  models  that  can  be  used  to  represent 
concrete,  namely,  material  type  5  (soil  and  crushable  foam),  material 
type  14  (soil  and  crushable  foam  failure),  material  type  16  (pseudo 
tensor),  material  type  25  (geological  cap  model),  material  type 
72RW3  (concrete  damage),  material  type  84  (winfrith  concrete),  and 
material  type  96  (brittle  damage).  Material  type  72RW3  (MAT_- 
CONCRETE_DAMAGE),  was  the  third  release  of  Karagozian  and  Case 
(K  and  C)  concrete  model.  It  is  a  plasticity-based  model,  using  three 
shear  failure  surfaces  and  including  damage  and  strain  rate  effects 
[8].  The  model  has  a  default  parameter  generation  function  based  on 
the  unconflned  compressive  strength  of  the  concrete  and  provides 
a  robust  representation  of  complex  concrete  laboratory  response  [9]. 
In  this  model,  the  stress  tensor  is  expressed  as  the  sum  of  the 
hydrostatic  stress  tensor  and  the  deviatoric  stress  tensor.  The 
hydrostatic  tensor  changes  the  concrete  volume  and  the  deviatoric 
stress  tensor  controls  the  shape  deformation. 

For  the  hydrostatic  stress  tensor,  the  compaction  model  is 
a  multi-linear  approximation  in  internal  energy.  Pressure  is  defined 
by 

p  =  C{s^)  +  yT{s^)E  (1) 


Fig.  2.  Finite  element  model  of  reinforced  concrete  columns. 


Fig.  3.  Pressure  versus  volumetric  strain  curve. 


Ao  >  0 


Fig.  4.  Strength  model  for  concrete  [8]:  (a)  failure  surfaces  in  concrete  material  model;  (b)  concrete  constitutive  model. 


where  E  is  the  internal  energy  per  initial  volume,  y  is  the  ratio  of 
specific  heats.  The  volumetric  strain,  (^v),  is  given  by  the  natural 
logarithm  of  the  relative  volume.  As  shown  in  Fig.  3  the  model 
contains  an  elastic  path  from  the  hydrostatic  tension  cut-off  to  the 
point  T  of  elastic  limit.  When  tension  stress  is  greater  than  the 
hydrostatic  tension  cut-off,  tension  failure  occurs.  When  the  volu¬ 
metric  strain  exceeds  the  elastic  limit,  compaction  occurs  and  the 
concrete  turns  into  a  granular  kind  of  material.  The  bulk  unloading 
modulus  is  a  function  of  volumetric  strain.  Unloading  occurs  along 
the  unloading  bulk  modulus  to  the  pressure  cut-off.  Reloading 
always  follows  the  unloading  path  to  the  point  where  unloading 
began,  and  continues  on  the  loading  path. 

A  three-curve  model  is  used  to  analyze  the  deviatoric  stress 
tensor,  as  shown  in  Fig.  4,  where  the  upper  curve  represents  the 
maximum  strength  curve,  the  middle  curve  is  the  initial  yield 
strength  curve  and  the  lower  curve  is  the  failed  material  residual 
strength  curve. 

In  order  to  consider  the  fact  that  under  higher  loading  rates 
concrete  exhibited  increased  strength,  a  dynamic  increase  factor 
(DIF),  the  ratio  of  the  dynamic  to  static  strength,  is  employed  in  this 
analysis.  The  expressions  proposed  by  Malvar  and  Crawford  [10,11] 
are  utilized.  The  DIF  for  the  concrete  compressive  strength  is  given  as: 
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Strain  Rate  (1/s) 


DIF 


f<30s-i 

£>30S“'‘ 


(2) 


where  e  is  the  strain  rate  in  the  range  of  30  x  10“®  to  300  s“^ 
£s  =  30  X  10“®  s“^  (static  strain  rate);  log  7s  =  6.156a:s-2; 
as  =  1/(5  +  9/c//co);  /co^lOMPa;  fc  is  the  static  compressive 
strength  of  concrete.  A  plot  of  the  formulae  employed  in  this 
study  for  the  DIF  of  concrete  in  compression  is  shown  in  Fig.  5(a). 

The  DIF  for  concrete  in  tension  is  given  by: 


DIF 


(e/fs)^  e<1.0s“^ 

P^e/esY^^  £  >  1.0  S“^ 


(3) 


where  e  is  the  strain  rate  in  the  range  of  10  ®s  ^  to  160 s 
=  10“®  s“^  (static  Strain  rate);  log =  65-2]  5  =  1/(1  +8/c//co); 
fco  =  10  MPa;/c  is  the  static  compressive  strength  of  concrete.  A  plot  of 
the  formulae  employed  in  this  study  for  the  DIF  of  concrete  in 
compression  is  shown  in  Fig.  5(b). 

Based  on  Fig.  5(a)  and  (b),  it  can  be  seen  that  the  tensile  response 
is  more  sensitive  to  strain  rate  than  compressive  response.  Therefore, 
different  rate  enhancements  are  included  in  “tension  and 
compression”  in  the  concrete  material  model  employed  in  this  study. 
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Fig.  5.  (a)  DIF  for  concrete  in  compression  [10].  (b)  DIF  for  concrete  in  tension  [10]. 
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In  the  numerical  model,  the  strain  rate  effect  is  incorporated  as 
follows.  At  any  given  pressure,  the  failure  surfaces  are  expanded  by 
a  rate  enhancement  factor  which  depends  on  the  effective  devia- 
toric  strain  rate,  as  shown  in  Fig.  6.  Let  rf  be  the  strain  rate 
enhancement  factor  and  p  the  pressure;  an  “unenhanced”  pressure 
p/rf  is  first  obtained,  then  the  unenhanced  strength  d(7(p/rf)  is 
calculated  for  the  specified  failure  surface.  Finally,  the  enhanced 
strength  is  given  by: 

A(Te  =  rfAff(p/rf)  (4) 


Strain  Rate  (1/s) 


Fig.  7.  DIF  for  reinforcement  [10]. 

Eqs.  (5)  and  (6)  are  valid  for  reinforcement  with  yield  stress 
between  290  and  710  MPa  and  for  strain  rates  between  10“^  s“^ 
and  225  s“\ 


2.4.  Application  of  loads  and  analysis  procedure 

The  behaviour  of  columns  subjected  to  blast  conditions  will  be 
influenced  by  the  fact  that  in  most  cases,  the  columns  would  have 
already  been  subjected  to  their  respective  service  prior  to  being 
exposed  to  blast  effects.  Therefore,  in  the  first  loading  stage,  gravity 
load  is  applied  via  slow  ramps  to  the  column,  while  in  the  second 
loading  stage,  blast  and  gravity  loads  are  applied  simultaneously. 
The  gravity  load  includes  both  the  dead  and  service  loads  acting  on 
the  column.  The  gravity  loads  are  assumed  to  be  in  the  range  of 
0.1— 0.4/^Ag  to  investigate  the  effects  of  gravity  load  on  the 
dynamic  response  and  residual  axial  capacity  of  reinforced  concrete 
columns  under  short  standoff  blast  conditions.  Eventually,  in  the 
post-blast  loading  stage,  axial  load  is  increased  until  the  column  is 
crushed,  as  shown  in  Fig.  8. 


Strength  is  equally  enhanced  along  any  radial  stress  path, 
including  uniaxial,  biaxial  and  triaxial  tension,  and  uniaxial  and 
biaxial  compression.  The  effective  strain  rate  versus  deviatoric 
strength  enhancement  is  given  by  a  LS-DYNA  define  curve  keyword. 


3.  Verification  of  finite  element  model 

Verification  of  the  finite  element  models  as  outlined  in  the 
above  procedure  is  carried  out  by  performing  the  analysis  of  several 
correlated  studies. 


2.3.2.  Reinforcement 

Steel  is  modelled  as  a  strain  rate  sensitive  uniaxial  elasto-plastic 
material  to  account  for  its  strain  rate  sensitivity  and  stress-strain 
history  dependence.  The  stress-strain  curve  is  assumed  to  be 
bilinear,  representing  an  elasto-plastic  behaviour  with  linear 
isotropic  hardening.  For  the  strain  rate  sensitivity,  the  expressions 
proposed  by  Malvar  and  Crawford  [10]  are  utilized.  A  plot  of  the 
proposed  formulae  is  shown  in  Fig.  7. 

The  yield  stress  of  reinforcement  is  represented  by: 

DIF  =  (5) 

where  a  =  a:fy,a:fy  =  0. 074-0. 04/y/414  and /y  is  the  static  yield 
strength  of  reinforcement  in  MPa. 

The  ultimate  stress  of  reinforcement  is  represented  by: 

DIF  =  (6) 

where  a  =  =  0. 019-0. 009/u/414,  and  /u  is  the  static 

ultimate  strength  of  reinforcement  in  MPa. 


3.1.  Dynamic  testing  on  concentrically  loaded  reinforced  concrete 
columns  confined  by  square  hoops 

The  transverse  reinforcement  is  modelled  using  beam  elements 
to  take  into  account  the  confinement  effect.  In  addition,  the 
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Fig.  8.  Loading  procedures  for  finite  element  analysis. 


concrete  material  model  has  one  parameter  (hi)  that  governs 
softening  in  compression.  The  axial  load-strain  curve  of  a  small 
reinforced  concrete  column  for  a  uniaxial  compression  test  con¬ 
ducted  by  Mander  et  al  was  reproduced  [12]  for  model  calibration. 

During  the  experiment,  the  column  was  loaded  concentrically  in 
a  DARTEC  10  MN  servo  hydraulically  controlled  testing  machine. 
Because  of  the  high  oil-pumping  capacity,  a  high  axial  strain  rate 
could  be  achieved.  For  the  selected  specimen,  the  strain  rate  was 
Sc  =  0.0167/s.  The  axial  strain  was  recorded  over  the  central 
400  mm  gauge  length  of  column  using  four  linear  potentiometers. 
The  axial  strain  plotted  is  defined  as  the  average  strain  occurring  in 
those  gauge  lengths  within  the  crushing  region  rather  than  the 
average  strain  of  all  gauge  lengths.  Therefore,  the  critical  axial 
strain  at  the  central  portion  of  the  column  from  the  numerical 
analysis  is  used  to  match  the  experimental  data.  The  experimental 
and  analytical  results  are  plot  in  Fig.  9.  Agreement  is  close  with  the 
predicted  peak  axial  load  exceeding  the  experimental  results  by 
around  8%. 

3.2.  Reinforced  concrete  columns  subjected  to  simulated 
blast  loading 


Specific  Impulse  (psi-sec) 

Fig.  10.  Comparison  of  numerical  and  experimental  residual  deflections. 


The  first  explosive  loading  laboratory  testing  program  at  the 
University  of  California,  San  Diego  utilizes  a  hydraulic-based  blast 
simulator  to  simulate  explosive  events  without  using  explosive 
materials  [13].  Several  tests  have  been  performed  to  investigate  the 
dynamic  response  of  the  reinforced  concrete  columns  when  sub¬ 
jected  to  impulsive  loads. 

The  dynamic  responses  of  the  test  specimens  subjected  to 
impulse  loads  of  0.77  psi-s,  1.76  psi-s,  1.9  psi-s  and  2.3  psi-s  were 
analyzed  using  the  proposed  finite  element  models.  Three  cases  of 
positive  duration  representative  of  the  typical  energy  dissipation 
time  of  a  close-in  explosion,  3  ms,  4  ms,  and  5  ms,  and  their 
respective  peak  pressure  were  utilized  in  the  analysis,  as  the 
detailed  peak  pressure  and  duration  for  the  corresponding  impulse 
loads  were  not  given. 

Fig.  10  shows  the  comparisons  of  residual  deformations  of 
numerical,  laboratory  and  field  test  results.  These  comparisons 
generally  show  a  good  agreement.  The  comparisons  show  that  the 
numerical  result  is  much  higher  than  the  laboratory  test  result 
when  subjected  to  an  impulse  load  of  1.76  psi-s.  While  for  the  other 
cases,  the  comparisons  generally  show  a  good  agreement.  Consid¬ 
ering  that  only  limited  data  are  available  and  the  unstable  character 
of  blast  test  results,  these  comparisons  are  considered  to  be 


reasonably  in  good  agreement.  Good  correlations  of  failure  mech¬ 
anisms  are  also  observed,  as  shown  in  Fig.  11.  From  both  field  test 
and  blast  simulator  test  results,  it  is  apparent  that  the  column  failed 
primarily  in  diagonal  shear  near  the  top  and  bottom  ends,  with  the 
central  portion  remaining  relatively  intact.  The  predicted  damage 
on  the  column  by  the  FE  model  is  shown  by  plotting  fringes  of 
effective  strain  which  is  used  for  measuring  the  overall  deformation 
at  one  point.  These  effective  strain  contours  reveal  the  strain 
localization  where  failure  propagates.  It  shows  that  the  failure  is 
localized  near  the  column  top  and  bottom  ends  due  to  diagonal 
shear  failure,  which  is  consistent  with  field  and  laboratory  test 
results. 

4.  Numerical  simulation  study 

4.1.  Numerical  simulation  matrix 

Using  the  finite  element  models  discussed  above,  numerical 
simulations  were  performed  to  evaluate  the  dynamic  response  in 
blast  situations,  and  to  further  estimate  the  residual  axial  capacity 
of  the  damaged  columns.  An  extensive  parametric  study  was 
carried  out  with  the  following  cases  considered  for  each 


Fig.  9.  Comparison  of  experimental  and  analytical  results. 
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Fig.  11.  Comparison  of  numerical  and  experimental  response  of  reinforced  concrete 
columns  subjected  to  impulsive  loads. 
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Fig.  12.  Column  geometry  and  reinforcement  details  of  the  simulation  matrix. 


parameter:  transverse  reinforcement  ratio  pv  =  0.12%  and  0.46%; 
longitudinal  reinforcement  pg  =  1.8%  and  3.2%;  column  aspect 
ratio  Ljb  =  9.8,  8  and  6;  axial  load  ratio  Pi/Zc^g  =  0-2,  0.3  and 

0.4.  The  parametric  study  was  carried  out  using  a  series  of  12 
columns,  labelled  as  series  A-L,  as  shown  in  Fig.  12.  The  series  B, 


D,  F,  FI,  J,  and  L  are  seismically  detailed  columns  for  which  the 
spacing  of  the  transverse  reinforcement  is  determined  in  accor¬ 
dance  with  the  requirement  in  the  ACl  318  code.  Table  1 
summarises  the  specimen  characteristics  of  the  simulation 
matrix. 
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Fig.  12.  (continued). 
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4.2.  Numerical  results  of  reinforced  concrete  columns  in  the 
dynamic  response  stage 

42 A.  Effect  of  transverse  reinforcement  ratio 

Blast  loadings  are  many  times  greater  than  conventional  loads. 
The  desired  ductile  flexural  behaviour  can  only  be  developed  when 
the  shear  capacity  exceeds  flexural  capacity.  The  ductility  capacity 
of  a  column  depends  on  the  amount  and  distribution  of  transverse 


reinforcement  within  the  plastic  hinge  region.  The  transverse 
reinforcement  increases  the  shear  capacity  of  the  column,  but  more 
importantly,  it  provides  conflnement  to  the  core  concrete  and 
lateral  restraint  against  buckling  of  the  longitudinal  reinforcement. 
Such  restraint  is  vital  for  reinforced  concrete  columns  which  have 
begun  to  crack  and  have  lost  the  majority  of  their  tensile  and 
flexural  capacity  but  still  need  to  bear  compressive  force  until 
ductile  failure  of  the  longitudinal  reinforcement  occurs.  Therefore, 


Table  1 

Specimen  characteristics  of  the  simulation  matrix. 


Column  type 

Cross  section 

Column  height 

Aspect  ratio 

Long,  steel 

Transverse  steel 

A 

355  X  355  mm  (14  X  Min) 

3480  mm  (137  in) 

9.8 

8T25  (pg  = 

3.2%) 

TlO-350  mm  (pv  =  0.12%) 

B 

355  X  355  mm  (14  X  Min) 

3480  mm  (137  in) 

9.8 

8T25  (pg  = 

3.2%) 

TlO-88  mm  (pv  =  0.46%) 

C 

355  X  355  mm  (14  X  Min) 

3480  mm  (137  in) 

9.8 

8T20  (pg  = 

1.8%) 

TlO-350  mm  (pv  =  0.12%) 

D 

355  X  355  mm  (14  X  Min) 

3480  mm  (137  in) 

9.8 

8T20  (pg  = 

1.8%) 

TlO-88  mm  (pv  =  0.46%) 

E 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T25  (pg  = 

3.2%) 

TlO-350  mm  (pv  =  0.12%) 

F 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T25  (pg  = 

3.2%) 

TlO-88  mm  (pv  =  0.46%) 

G 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T20  (pg  = 

1.8%) 

TlO-350  mm  (pv  =  0.12%) 

H 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T20  (pg  = 

1.8%) 

TlO-88  mm  (pv  =  0.46%) 

1 

355  X  355  mm  (14  X  Min) 

2130  mm  (84  in) 

6 

8T25  (pg  = 

3.2%) 

TlO-350  mm  (pv  =  0.12%) 

J 

355  X  355  mm  (14  X  Min) 

2130  mm  (84  in) 

6 

8T25  (pg  = 

3.2%) 

TlO-88  mm  (pv  =  0.46%) 

1< 
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2130  mm  (84  in) 

6 

8T20  (pg  = 

1.8%) 

TlO-350  mm  (pv  =  0.12%) 

L 

355  X  355  mm  (14  X  Min) 

2130  mm  (84  in) 

6 

8T20  (pg  = 

1.8%) 

TlO-88  mm  (pv  =  0.46%) 
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Fig.  13.  Effect  of  transverse  reinforcement  ratio  on  the  displacement  response  of  columns. 
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Fig.  14.  Effect  of  axial  load  ratio  on  the  displacement  response  of  columns  with  a  high  transverse  reinforcement  ratio. 


transverse  reinforcement  is  expected  to  have  a  significant  infiuence 
on  the  failure  mode  of  the  columns  under  blast  loading  and  their 
respective  blast  resistance.  Fig.  13  shows  the  effect  of  the  transverse 
reinforcement  ratio  on  the  displacement  response  of  the  columns. 
In  these  figures,  the  horizontal  axis  represents  the  magnitude  of  the 
impulse,  which  is  obtained  from  the  CONWEP  software;  the  vertical 
axial  yr/L  represents  the  displacement  to  height  ratio,  where  is 
the  residual  lateral  displacement  at  mid-height  and  L  is  the  clear 
height  of  the  column.  The  dotted  lines  denote  the  numerical  results, 
and  the  solid  lines  represent  the  fitted  trend  lines.  Series  A  and  C 
columns  are  conventional  columns  which  are  mainly  designed  for 
gravity  loads  and  are  weak  in  their  shear  capacity.  Series  B  and  D 
are  seismically  detailed  columns.  It  can  be  seen  that  under  the  same 
impulsive  loads,  the  shear-critical  columns  have  a  much  larger 
deflection  than  the  seismically  detailed  columns.  This  is  consistent 
with  the  predication  that  implementation  of  seismic  detailing  can 
significantly  reduce  the  degree  of  direct  damage  due  to  blast  loads 


and  consequently  improve  the  blast  resistance  of  the  columns.  The 
improvement  becomes  more  remarkable  in  cases  with  severe 
impulsive  loads. 

4.2.2.  Effect  of  axial  load 

Before  the  occurrence  of  a  blast  incident,  the  gravity  load  is 
already  imposed  on  the  column.  This  will  infiuence  the  behaviour 
of  the  column  under  blast  conditions.  Fig.  14  illustrates  the  effect  of 
axial  load  on  the  displacement  response  of  the  columns  with  a  high 
transverse  reinforcement  ratio.  As  the  figures  demonstrate,  when 
the  impulsive  loading  and  corresponding  deformation  are  small, 
the  mid-height  displacement  of  the  column  with  larger  axial  loads 
is  slightly  smaller.  This  is  due  to  the  fact  that  with  an  increase  in  the 
applied  axial  load  on  columns,  it  would  result  in  an  increase  in  its 
moment  capacity  and  its  nominal  shear  strength.  Flowever,  the 
decrease  in  mid-height  displacement  would  only  occur  before  the 
impulse  and  its  corresponding  displacement  reach  a  critical  value. 


Fig.  15.  Effect  of  axial  load  ratio  on  the  displacement  response  of  columns  with  a  low  transverse  reinforcement  ratio. 
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Fig.  16.  Effect  of  longitudinal  reinforcement  ratio  on  the  displacement  response  of  columns. 
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Fig.  17.  Effect  of  column  aspect  ratio  on  the  displacement  response  of  columns. 


Once  this  critical  value  is  exceeded,  the  mid-height  displacement 
would  increase  greatly  with  increasing  axial  load.  This  is  expected 
for  columns  with  flexural  behaviour.  When  columns  experience 
large  deflection  and  plastic  hinges  formation  occurs  at  mid-span 
and  flxed  ends,  axial  loads  will  amplify  the  lateral  deflection  and 
internal  moment  due  to  the  P-A  effect.  As  the  deflection  increases, 
the  column  will  transit  from  a  gradual  stiffness  and  strength 
degradation  to  a  rapid  loss  of  strength  due  to  the  buckling  of  the 
longitudinal  reinforcement.  This  explains  the  change  in  the  slope  of 
the  curve  from  gentle  to  steep  as  the  axial  loads  increase. 

Fig.  15  illustrates  the  effect  of  the  axial  load  on  the  displacement 
response  of  the  columns  with  a  low  transverse  reinforcement  ratio. 
When  impulsive  loading  and  the  corresponding  deformation  are 
small,  the  mid-height  displacement  of  the  column  with  larger  axial 
loads  is  slightly  smaller,  as  shown  in  the  flgure.  Columns  with 
higher  axial  loads  collapsed  when  the  impulse  loading  and  corre¬ 
sponding  deformation  were  small  due  to  the  brittle  behaviour  of 
shear-critical  columns.  Thus,  there  is  less  data  for  columns  with 
higher  axial  loads  shown  in  the  flgure.  It  is  believed  that  axial  load 
failure  occurred  immediately  after  the  shear  failure;  while  for  the 
columns  having  lower  axial  loads,  the  collapse  only  occurred  when 
the  impulsive  loading  and  corresponding  deformation  are  relatively 
large,  although  shear  failure  had  occurred  at  a  smaller  displacement 
ratio.  This  may  be  due  to  the  following  factor:  the  longitudinal 
reinforcement  will  support  a  portion  of  the  axial  load  up  to 
a  maximum  load  defined  by  either  the  buckling  (bulking  length 
equal  to  the  spacing  of  ties)  or  the  plastic  capacity  of  the  rein¬ 
forcement  bars;  when  the  axial  loads  were  less  than  this  capacity, 
the  axial  load  collapse  didn’t  occur  until  the  ductile  failure  of  the 
reinforcement.  Therefore,  the  columns  under  low  axial  loads  tend 
to  have  a  more  ductile  response  than  those  under  high  axial  loads. 


Similar  flndings  have  been  reported  for  both  reinforced  concrete 
columns  and  shear  walls  during  seismic  investigations  [14]. 

4.2.3.  Effect  of  longitudinal  reinforcement  ratio 

As  the  longitudinal  reinforcement  ratio  increases,  both  the 
ultimate  moment  capacity  and  axial  capacity  increase.  As  a  result, 
an  increase  in  longitudinal  reinforcement  ratio  would  improve  the 
blast  resistance  of  the  column  with  flexural  behaviour.  Fig.  16(a) 
illustrates  the  effect  of  the  longitudinal  reinforcement  ratio  on  the 


Fig.  19.  Effect  of  axial  load  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


displacement  response  of  the  columns  with  a  high  transverse 
reinforcement  ratio.  The  results  show  that  the  mid-height 
displacement  decreases  with  increasing  longitudinal  reinforce¬ 
ment  ratio  as  expected.  However,  an  increase  in  longitudinal  rein¬ 
forcement  may  not  always  improve  the  blast  resistance.  When  the 
flexural  capacity  of  the  column  exceeds  its  shear  capacity,  it  may 
lead  to  a  shift  from  a  ductile  flexural  failure  to  a  brittle  shear  failure 
mode.  In  such  cases,  the  mid-height  displacement  will  increase 
with  increasing  longitudinal  reinforcement  ratio.  Fig.  16(b)  shows 
the  effect  of  longitudinal  reinforcement  for  columns  with  a  low 
transverse  reinforcement  ratio.  In  this  flgure,  series  A  columns  are 
critical  in  shear  due  to  inadequate  transverse  reinforcement,  and 
series  C  columns  are  critical  in  flexure  due  to  the  fact  that  they  have 
lesser  flexural  capacity  than  shear  capacity  because  of  the  low 
longitudinal  reinforcement  ratio.  As  the  flgures  demonstrate,  when 
the  impulsive  loading  and  corresponding  deformation  are  small, 
the  mid-height  displacement  of  series  A  column  is  larger  than  that 
of  the  series  C  column.  This  is  consistent  with  the  prediction  above. 
However,  when  the  impulsive  loading  and  corresponding  defor¬ 
mation  are  large,  the  mid-height  displacement  of  series  C  column  is 


larger  than  that  of  series  A  column.  This  is  believed  to  result  from 
the  P-A  effect  for  flexural  columns  when  the  deflection  is  large. 
Flexural  failure  is  preferred  than  shear  failure  because  an  extended 
plastic  response  is  provided  prior  to  the  collapse  of  column. 
However,  it  is  noteworthy  that  the  deformation  in  flexure  beyond 
a  certain  limit  will  jeopardize  the  ability  of  the  column  to  carry 
vertical  loads  due  to  the  P-A  effect.  Therefore,  to  improve  the  blast 
resistance  and  assure  a  ductile  response,  sections  need  to  be 
designed  so  that  the  flexural  capacity  is  less  than  the  capacity  of 
non-ductile  failure  mechanism  and  the  deflection  must  be 
controlled  to  prevent  column  instability  due  to  the  P-A  effect. 

4.2.4.  Effect  of  column  aspect  ratio 

In  this  study,  the  cross  section  is  kept  constant  for  all  the  column 
specimens.  A  decrease  in  column  aspect  ratio  means  a  decrease  of 
overall  height.  As  a  result,  when  subjected  to  the  same  intensity  of  blast 
loads,  the  moment  and  shear  demand  will  decrease  with  decreasing 
column  aspect  ratio.  The  effect  of  aspect  ratio  on  the  displacement 
response  is  shown  in  Fig.  17.  It  is  observed  that  under  the  same 
impulsive  loading  the  displacement  to  height  ratio  decreases  with 


Fig.  20.  The  v—[{yr/L)  x  (P/fAg)]  curves  of  columns. 


Fig.  21.  Effect  of  longitudinal  reinforcement  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


a  reduced  aspect  ratio  as  expected.  Based  on  the  figures,  the  results  are 
similar  for  both  shear-critical  and  flexural  columns. 

4.3.  Numerical  results  of  residual  axial  capacity  of  the  blast 
damaged  reinforced  concrete 

In  the  post-blast  analysis  stage,  the  axial  load  is  gradually 
increased  by  applying  a  rigid  plate  attached  to  the  top  end  of  the 
column  using  the  displacement  mode  to  capture  both  the  residual 
axial  capacity  and  the  softening  portion  of  the  loading  curve.  One 
typical  curve  is  shown  in  Fig.  18.  It  is  noted  that  in  the  blast  loading 
stage,  due  to  the  inertia  effect,  the  axial  load  supported  by  the 
column  is  not  constant  but  fluctuates  along  the  deformation  of  the 
column.  The  sudden  increase  of  the  axial  force  at  a  displacement  of 
75  mm  could  be  due  to  the  steel  reinforcing  bars  reaching  their 
strain  hardening  stage.  The  increase  in  the  strength  of  the  steel 
reinforcement  at  this  point  in  time  would  result  in  an  increase  in 
the  columns  axial  load  carrying  capacity.  The  ultimate  state  of 
reinforced  concrete  columns  has  often  been  defined  by  some 


researchers  [15]  as  a  state  of  vanishing  axial  capacity  to  sustain  the 
dead  and  live  loads  (long-term  load),  which  is  indicated  as  the  point 
of  axial  failure  (Fig.  19). 

4.3.1.  Post-blast  damage  evaluation 

In  this  study,  the  residual  axial  capacity  level  of  a  blast  damaged 
column  was  evaluated  by  the  ratio  of  residual  axial  strength,  which 
was  defined  by  the  following  equation: 

=  (Pr-PO/iPmax-Pl)  (7) 

where  Pmax  is  the  axial  capacity  of  the  undamaged  columns,  Pr  is 
the  residual  axial  capacity  of  the  damaged  columns  and  Pl  is  the 
long-term  axial  load. 

When  the  column  is  undamaged,  Pr  =  Pmax,  the  value  of  v  is  1 ; 
when  the  column  has  lost  the  ability  to  sustain  the  long-term  axial 
load,  Pr  =  Pl,  the  value  of  v  just  reaches  zero,  referring  as  the 
ultimate  state  of  the  column.  As  for  the  performance  indicator,  the 
previously  defined  displacement  to  height  ratio  (yr/I)  is  used. 


Fig.  22.  Effect  of  transverse  reinforcement  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


Fig.  23.  Effect  of  column  aspect  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


4.3.2.  Effect  of  axial  load  ratio 

Fig.  19  shows  the  effect  of  axial  load  ratio  on  the  residual  axial 
capacity  of  the  columns  at  various  degrees  of  deformation  level. 
The  results  show  that  at  the  same  mid-height  displacement  ratio, 
the  ratio  of  residual  axial  capacity  is  smaller  in  the  case  of  larger 


axial  loads.  As  illustrated  by  the  figure,  for  columns  with  lower  axial 
loads,  the  axial  load  failure  tends  to  occur  at  a  relatively  large 
displacement  ratio;  for  columns  with  a  larger  axial  load,  axial  load 
failure  tends  to  occur  at  a  smaller  displacement  ratio.  This  indicates 
that  the  displacement  ratio  at  axial  load  failure  is  inversely  related 


Fig.  24.  Comparison  of  numerical  results  with  the  proposed  curves. 


to  the  magnitude  of  axial  load.  It  is  also  noted  that  for  columns  with 
a  low  transverse  reinforcement  ratio,  the  effect  of  axial  load  is  more 
critical.  This  is  believed  to  be  due  to  the  significance  of  confinement 
during  axial  compression. 

Based  on  the  observation,  a  new  term  (yr/T)  x  (Pi/f^Ag)  was 
introduced.  Fig.  20  plots  the  relation  between  v  and 
(yr/L)  X  (Pi/f^Ag)  for  the  columns.  It  can  be  seen  that  this  new  term 
is  capable  of  reflecting  the  influence  of  the  axial  load  ratio  on  the 
ratio  of  residual  axial  capacity.  It  also  highlights  the  effects  of  other 
parameters  and  enables  them  to  be  studied  based  on  the  compar¬ 
ison  of  V  versus  (yr/L)  x  (Pi/f/Ag)  curves. 

4.3.3.  Effect  of  longitudinal  reinforcement  ratio 

It  is  usually  assumed  that  longitudinal  reinforcement  will 
support  a  portion  of  the  axial  load  up  to  a  maximum  load  defined 
by  either  the  buckling  or  the  plastic  capacity  of  the  reinforcing  bars. 
In  most  cases  the  column  collapse  is  related  to  the  increase  of  axial 
load  carried  by  the  longitudinal  reinforcing  bars  and  their  deteri¬ 
oration  of  compressive  strength.  Fig.  21  shows  the  effect  of  the 
longitudinal  reinforcement  ratio  on  the  residual  axial  capacity  ratio 
of  the  columns.  The  results  indicate  that  the  ratio  of  residual  axial 
capacity  is  generally  larger  when  the  longitudinal  reinforcement 
ratio  increases. 

4.3.4.  Effect  of  transverse  reinforcement  ratio 

Fig.  22  shows  the  effect  of  the  transverse  reinforcement  ratio  on 
the  residual  axial  capacity  of  the  columns  at  various  degrees  of 
deformation  level.  The  results  show  that  the  ratio  of  residual  axial 
capacity  of  columns  with  a  low  transverse  reinforcement  ratio  is 
significantly  less  than  that  of  the  column  with  a  high  transverse 
reinforcement  ratio.  From  these  figures,  it  is  observed  that  under 
similar  axial  loading  conditions,  axial  load  failure  tends  to  occur  at 
a  relatively  large  displacement  ratio  for  columns  with  high  trans¬ 
verse  reinforcement  ratio.  In  contrast,  columns  with  a  low  trans¬ 
verse  reinforcement  ratio  will  tend  to  fail  under  axial  loads  at 
a  smaller  displacement  ratio.  This  suggests  that  the  displacement 
ratio  at  axial  load  failure  is  directly  related  to  the  transverse  rein¬ 
forcement  ratio. 

4.3.5.  Effect  of  column  aspect  ratio 

Fig.  23  illustrates  the  effect  of  the  column  aspect  ratio  on  the 
residual  axial  strength  ratio  of  the  columns.  The  results  show  that  at 
a  high  transverse  reinforcement  ratio,  the  residual  axial  capacity 
ratio  increases  with  the  reduction  in  aspect  ratio.  For  columns  with 
a  low  transverse  reinforcement  ratio,  this  effect  is  not  very  clear  due 
to  the  scatter  in  the  results. 

4.4.  Proposed  formulae  for  determining  the  residual  axial  capacity 
ratio  of  blast  damaged  reinforced  concrete  columns 

The  parametric  study  carried  out  revealed  the  significance  of  the 
parameters  that  affect  the  residual  axial  strength  of  the  blast 
damaged  reinforced  concrete  column.  A  formula  was  derived 
through  multivariable  regression  analysis  in  terms  of  various 
parameters  to  predict  the  residual  axial  capacity  ratio  based  on  the 
mid-height  displacement  to  height  ratio  by  fitting  the  parametric 
study  results  and  is  as  follows: 

V  =  [73.65pv  +  8.465pg-0.020879(L/fa)+0.104] 

g[89284.22p^-1308.64221pg-9.684203(L/b)-382.12](yr/L)(PL//cA)  (g) 

A  few  examples  presenting  the  comparison  of  the  proposed 
equation  with  the  analytical  results  are  shown  in  Fig.  24.  The  solid 
line  in  each  plot  denotes  the  proposed  equation,  while  the  dotted 


line  represents  the  numerical  analysis  results.  It  is  observed  that 
for  most  parts,  the  proposed  curves  are  close  to  the  analytical 
results. 

5.  Summary  and  conclusions 

Based  on  the  results  of  the  parametric  study,  the  following  main 
conclusions  can  be  drawn. 

-  The  numerical  results  show  that  the  use  of  seismic  detailing 
techniques  can  significantly  reduce  the  degree  of  direct  blast- 
induced  damage  and  subsequent  collapse  of  the  reinforced 
concrete  columns. 

-  Comparisons  of  the  deterioration  of  the  axial  strength  under 
different  axial  load  ratios  indicate  that  the  ratio  of  residual  axial 
strength  is  smaller  under  larger  long-term  axial  load.  The  effect 
of  axial  load  ratio  is  more  critical  in  the  case  of  columns  with 
a  low  transverse  reinforcement  ratio. 

-  The  results  indicate  that  the  ratio  of  residual  axial  capacity 
generally  increases  with  an  increase  in  the  longitudinal  rein¬ 
forcement  ratio. 

-  The  numerical  results  indicate  that  the  residual  axial  capacity 
ratio  increases  with  a  reduction  in  the  aspect  ratio  for  columns 
with  a  high  transverse  reinforcement  ratio. 

-  A  formula  was  derived  by  fitting  the  parametric  study  results, 
in  terms  of  various  parameters  to  predict  the  residual  axial 
capacity  ratio  based  on  the  mid-height  displacement  to  height 
ratio.  The  comparison  of  the  proposed  equation  with  the 
analytical  results  shows  that  the  proposed  curve  well  repre¬ 
sents  the  tendency  with  the  variation  of  the  parameters.  Future 
experimental  investigation  of  residual  axial  strength  of  the 
blast  damaged  reinforced  concrete  columns  is  needed  to 
further  supplement  the  limited  data  set  used  to  develop  the 
proposed  equation. 
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Fig.  4.  Strength  model  for  concrete  [S]:  {a)  failure  surfaces  in  concrete  material  model;  (b)  concrete  constitutive  model. 


3.  Verification  of  finite  element  model 
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Fig.  8.  Loading  procedures  for  finite  element  analysis. 
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The  first  explosive  loading  laboratory  testing  program  at  the 
University  of  California,  San  Diego  utilizes  a  hydraulic- based  blast 
simulator  to  simulate  explosive  events  without  using  explosive 
materials  [13].  Several  tests  have  been  performed  to  investigate  the 
dynamic  response  of  the  reinforced  concrete  columns  when  sub¬ 
jected  to  impulsive  loads. 

The  dynamic  responses  of  the  test  specimens  subjected  to 
impulse  loads  of  0.77  psi-s,  1.76  psi-s,  1.9  psi-s  and  2.3  psi-s  were 
analyzed  using  the  proposed  finite  element  models.  Three  cases  of 
positive  duration  representative  of  the  typical  energy  dissipation 
time  of  a  close-in  explosion,  3  ms,  4  ms,  and  5  ms,  and  their 
respective  peak  pressure  were  utilized  in  the  analysis,  as  the 
detailed  peak  pressure  and  duration  for  the  corresponding  impulse 
loads  were  not  given. 
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Fig.  9.  Compiriiion  of  o;tporimjL-nLal  ind  analylitil  noifulLii. 
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Fig.  10,  Comparison  of  numerical  and  experimental  residual  deflections. 
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Column  geometry  and  reinforcement  details  of  the  simulation  matrix. 


Table  1 

Specimen  characteristics  of  the  simulation  matrix. 


Column  type 

Cross  section 

Column  height 

Aspect  ratio 

Long,  steel 

Transverse  steel 

A 

355  x355  mm  (14x14  in) 

3480  mm  (137  in) 

9.8 

8T25  (pg  =  3.2%) 

T1 0-350  mm(pv  =  0.12%) 

B 

355  X  355  mm  (14  x  14  in) 

3480  mm  (137  in) 

9.8 

8T25  (pg  =  3.2%) 

TlO-88  mm  (pv  =  0.46%) 

C 

355  X  355  mm  (14  x  14  in) 

3480  mm  (137  in) 

9.8 

8T20  (pg=h8%) 

T1 0-350  mm(pv  =  0.12%) 

D 

355  X  355  mm  (14  x  14  in) 

3480  mm  (137  in) 

9.8 

8T20  (pg=1.8%) 

T1 0-88  mm  (pv  =  0.46%) 

E 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T25  (pg  =  3.2%) 

T1 0-350  mm(pv  =  0.12%) 

F 

355  X  355  mm  (14  X  Min) 

2840  mm  (112  in) 

8 

8T25  (pg  =  3.2%) 

TlO-88  mm  (pv  =  0.46%) 

G 

355  X  355  mm  (14  x  14  in) 

2840  mm  (112  in) 

8 

8T20  (pg=1.8%) 

T1 0-350  mm(pv  =  0.12%) 

H 

355  X  355  mm  (14  x  14  in) 

2840  mm  (112  in) 

8 

8T20  (pg=1.8%) 

T1 0-88  mm  (pv  =  0.46%) 

I 

355  X  355  mm  (14  X  Min) 

2130  mm  (84  in) 

6 

8T25  (pg  =  3.2%) 

T1 0-350  mm(pv  =  0.12%) 

J 

355  X  355  mm  (14  X  Min) 

2130  mm  (84  in) 

6 

8T25  (pg  =  3.2%) 

TlO-88  mm  (pv  =  0.46%) 

I< 

355  X  355  mm  (14  x  14  in) 

2130  mm  (84  in) 

6 

8T20  (pg=1.8%) 

T1 0-350  mm(pv  =  0.12%) 

L 

355  X  355  mm  (14  x  Min) 

2130  mm  (84  in) 

6 

8T20  (pg=1.8%) 

T1 0-88  mm  (pv  =  0.46%) 
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Effect  of  transverse  reinforcement  ratio  on  the  displacement  response  of  columns. 


Effect  of  axial  load  ratio  on  the  displacement  response  of  columns  with  a  high  transverse  reinforcement  ratio. 
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Fig.  15.  Effect  of  axial  load  ratio  on  the  displacement  response  of  columns  with  a  low  transverse  reinforcement  ratio. 
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Fig.  16.  Effect  of  longitudinal  reinforcement  ratio  on  the  displacement  response  of  columns. 
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Fig.  17.  Effecl  of  column  aspect  ratio  on  the  displacement  response  of  columns. 
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4.3.1.  Post-blast  damage  evaluation 

In  this  study,  the  residual  axial  capacity  level  of  a  blast  damaged 
column  was  evaluated  by  the  ratio  of  residual  axial  strength,  which 
was  defined  by  the  following  equation: 

=  (Pr-PL)/(Pmax-PL)  (7) 

where  P^ax  is  the  axial  capacity  of  the  undamaged  columns,  Pj-  is 
the  residual  axial  capacity  of  the  damaged  columns  and  Pl  is  the 
long-term  axial  load. 

When  the  column  is  undamaged,  Pj-  =  Pmax,  the  value  of  v  is  1 ; 
when  the  column  has  lost  the  ability  to  sustain  the  long-term  axial 
load,  P,-  =  Pl,  the  value  of  v  just  reaches  zero,  referring  as  the 
ultimate  state  of  the  column.  As  for  the  performance  indicator,  the 
previously  defined  displacement  to  height  ratio  (yr/b)  is  used. 
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Fig.  19.  Effecl  of  axial  load  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 
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Fig.  20.  The  v—[{y^/V)  x  (P//cJ^g)]  curves  of  columns. 


Fig,  21.  EITecl  of  longitudinal  reinforcement  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


Fig.  22.  Effect  of  transverse  reinforcement  ratio  on  the 


of  residual  axial  capacity  of  the  blast  damaged  columns. 
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Fig.  23.  Effect  of  column  aspect  ratio  on  the  ratio  of  residual  axial  capacity  of  the  blast  damaged  columns. 


The  parametric  study  carried  out  revealed  the  significance  of  the 
parameters  that  affect  the  residual  axial  strength  of  the  blast 
damaged  reinforced  concrete  column.  A  formula  was  derived 
through  multivariable  regression  analysis  in  terms  of  various 
parameters  to  predict  the  residual  axial  capacity  ratio  based  on  the 
mid-height  displacement  to  height  ratio  by  fitting  the  parametric 
study  results  and  is  as  follows: 

V  =  [73.65pv  +  8.465pg-0.020879(L/£i)  +  0.104] 

p[89284.22(Ov-1308.6422t(Og-9.684203(t/fe)-382.12](yr/t)(PL/fA)  ("gl 
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.  Comparison  of  numerical  results  with  the  proposed  curves. 


NANYANG 

TECHNOLOGICAL 

UNIVERSITY 


(y,/L)‘(PL/f'cAg) 


5.  Summai^  and  conclusions 


Based  on  the  results  of  the  parametric  study,  the  following  main 
conclusions  can  be  drawn. 


-  The  numerical  results  show  that  the  use  of  seismic  detailing 
techniques  can  significantly  reduce  the  degree  of  direct  blast- 
induced  damage  and  subsequent  collapse  of  the  reinforced 
concrete  columns. 

-  Comparisons  of  the  deterioration  of  the  axial  strength  under 
different  axial  load  ratios  indicate  that  the  ratio  of  residual  axial 
strength  is  smaller  under  larger  long-term  axial  load.  The  effect 
of  axial  load  ratio  is  more  critical  in  the  case  of  columns  with 
a  low  transverse  reinforcement  ratio. 

-  The  results  indicate  that  the  ratio  of  residual  axial  capacity 
generally  increases  with  an  increase  in  the  longitudinal  rein¬ 
forcement  ratio. 
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5.  Summary  and  conclusions 


Based  on  the  results  of  the  parametric  study,  the  following  main 
conclusions  can  be  drawn. 

-  The  numerical  results  indicate  that  the  residual  axial  capacity 
ratio  increases  with  a  reduction  in  the  aspect  ratio  for  columns 
with  a  high  transverse  reinforcement  ratio. 

-  A  formula  was  derived  by  fitting  the  parametric  study  results, 
in  terms  of  various  parameters  to  predict  the  residual  axial 
capacity  ratio  based  on  the  mid-height  displacement  to  height 
ratio.  The  comparison  of  the  proposed  equation  with  the 
analytical  results  shows  that  the  proposed  curve  well  repre¬ 
sents  the  tendency  with  the  variation  of  the  parameters.  Future 
experimental  investigation  of  residual  axial  strength  of  the 
blast  damaged  reinforced  concrete  columns  is  needed  to 
further  supplement  the  limited  data  set  used  to  develop  the 
proposed  equation. 
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